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Abstract The irreducible tensor operators and their tensor products employing Racah al- 
gebra are studied. Transformation procedure of the coordinate system operators act on are 
introduced. The rotation matrices and their parametrization by the spherical coordinates of 
vector in the fixed and rotated coordinate systems are determined. A new way of calculation 
of the irreducible coupled tensor product matrix elements is suggested. As an example, the 
proposed technique is applied for the matrix element construction for two electrons in a field 
of a fixed nucleus. 
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1 Introduction 

in 



The main aim of present work is to parametrize irreducible matrix representation of either 
SO(3) or SU (2) group by the coordinates of S 2 x S 2 , where S 2 denotes the unit 2-dimension 
sphere. The motivation is grounded on the following occasions: (i) the difficulties in theoret- 
ical atomic spectroscopy arising through multiple integrals of Af-electron angular parts; (ii) 
inconvenience of Wigner-Eckart theorem application for irreducible tensor operator matrix 
elements on the basis of functions, expressed in terms of Wigner D-function. 
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In theoretical atomic physics the algorithms of matrix element calculation for atomic 
quantities on the basis of many-electron wave functions are well known and widely used 
Q]|2). The construction of matrix element is based on the structure of many-electron func- 
tion, which is represented by coupled tensor product of one-electron eigenstates. The latter 
formulation leads to complicated TV-electron angular parts and various techniques, in order 
to simplify the calculation of many-electron matrix element (3EU5). According to Racah 
[6 7 1, the basis eigenstates <j)^ behave in the same way as the spherical tensor operators 
T*; by means of transformations under irreducible matrix representations. The eigenstates 
0,^ of central-field atomic Hamiltonian are usually enunciated by spherical functions on 
SO(3)/SO(2) quotient group, multiplied by 2 x 1 spin matrices. However, in [8] it was 
showed that ^fe) = (j) (Q^,(p^) can be expressed by Wigner D (X2 5 ) functions on 517(2), 
where Q.^ = (3>e , ©£ , 0) = {(p^ + n/2, , 0) . Furthermore, Bhatia et al. in (9) constructed 
two-electron wave function <P{x\,X2) in terms of the spherical functions D(Q), where Q, 
as usually, denotes Euler angles (<P, ©, f), i.e., rotation on 5 2 from x\ to X2 and vice verse. 
Unfortunately, it is evident that for this basis the Wigner-Eckart theorem can not be applied 
directly, what leads to necessity to expand D (Q ) over x\ ,X2- 

In this work we start from notations in |T), [10|. We express rotation matrix by the 
spherical coordinates of vector in the fixed and rotated coordinate systems, rather than the 
rotation angles in an explicit form (Sec.[3]and Sec.[4]l. Following this route, we demonstrate 
technique of matrix element construction, when 2AMntegral (over the spherical coordinates 
0* , <p^ with E, = 1, 2, . . . ,N) is reduced up to a double one (Sec. |SJ, The technique is based 
on studied properties of integrity for obtained spherical functions (Sec. [5} and proposed 
transformation coefficients, called the rotated Clebsch-Gordan coefficients (CGC) or simply 
RCGC(Sec.l6Tl>. 



2 Preliminaries 

The well known transformation formula for &-rank spherical tensor operators T^,qE [—k, +k] , 
is given by 

T*(K 2 ) = 11%, (1) 

where the symbols K\ and Ki show that the corresponding tensor operator is defined in the 
fixed and rotated coordinate system, respectively. The generalized spherical function Tr , 
on 50(3) or SU{2) is of the form ED 



D k qq , (Q) = a (k,q,q') e i (^+«" p ) {cos Q©)} 2 ^ (k,q,q') {tan (|0)} 2 ^«+«' , (2) 



a (k, q,q')= i« -1 ^(k + q)l(k-q)\(k + q')\(k-q% (3) 

(-l) p 

b p {k,q,q')= \ — , (4) 

p[(p + q'-q)\(k + q-p)\(k-q'-p)\ 

where 4>, f 6 [0, 2k] and ® 6 [0, 7t\. In this work the standard phase system is used in which 
complex and Hermitian conjugate operators accordingly denote D* , = (— \) q ~ q ' D k _ q and 
= (—l) k ^ q T k _ q . It follows from Eq. (TJl, the rotated coordinate system K2 depends on the 
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fixed system ^1 and the rotation angles £2. Such a dependency can be reformulated in a 
different way: the rotation angles Q depend on the coordinates of given (or known) point, 
located on S 2 in the fixed and rotated coordinate systems. These spherical coordinates will 
be denoted x\ (in K\) and X2 (in K2). On the other hand, the spherical function D, which 
depends on Q, can also be expressed as a function of the variables x\, xq,. In the next section 
the latter dependency in an explicit form is presented. 



3 The geometry of rotation angles 

In order to express the generalized spherical function D by the coordinates x\ , x.2 of vector 
in the fixed and rotated coordinate system, the following geometry is defined. Suppose we 
have a map Q : S 2 x S 2 h-> SO (3). Its representation onto the linear space of vectors r, = 
ri/\ri\ = (xi yi n ) T € R 3 (|r,-| = const. Mi £ Z + - positive integers) is associated in the 
following way 

r 2 =R(n)r u (5) 
Xi = sin 8, cos 9/ , y,- = sin 0, sin 9,- , Zi = cos 0, , 

where the coordinate system K\ is rotated (the ZXZ convention) by the rotation matrix R of 
SO (3) I10I11I . Then the system of three equations in Eq. (f5]l is rewritten in the form 

(X2 = —x 1 ] sin 9 — u! cos 9 , 
y2 = x\ cos 9 — u' sin 9 , (6) 
Z2 =y' 1 sin0+z' 1 cos0, 

where 

m' = y'j cos 9 — z'x sin 9 , (7) 

!x' x = x\ cosf — yi sinf, 
y'j =x 1 sinf / +y 1 cosf , (8) 
z\=zi. 

The parameter f is chosen optionally in the range [0, 2k] . Partial solutions of subsystems 
(x2,Z2) and (y2,Z2) for 9,(p are substituted then in Eq. QJl. Optimal values of f are found 
by solving variational equation, varying obtained spherical function with respect to f (for 
details see Appendix lAll. We attain, that solutions in K satisfy equation 

sin 0i cos (<?>! +«P) =0. (9) 

The solutions are: 

1. cos (91 +f) = ¥ = -91 + a'f + CT7r + 27r«', n' e Z+, c' = ±1, a = 0, ±1. 

2. sin0i =0 Vfe [0,2w]. 

First of all let us analyze item (TTJ - the situation when sin 9\ ^ 0. Then y\ = a' sin 9\ , if 
(7 = and y'j = — CT'sin0i if d 1 = ±1 (see Eq. ©). Secondly, when item l[2} is valid, we 
obtain that 9\ = or 0j = K. The parameter f then could be of any value in the range 
[0, 2k] . It is clear that the solution sin 9\ = (item (O) is a particular case of item if the 
angle f is chosen to be equal to f = — 91 +0"'^ +CT^ + 2^«'. Since in item lf2]i the solution 
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Table 1 The values for the parameters a, j3, % 8, n 
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f can have arbitrary values in the range [0, 2k], we choose it equal to the solution given by 
item Q. Consequently, the solutions of the system in Eq. © are 

4> = <p 2 + ce|, = jS (0i - ye 2 )+2Kn, f = -91 + s| +2^n', (10) 

where ra' 6 Z + and the values for a, jS, y, 5, n are presented in Table [T] The function 
il{xi,X2) may be expanded into several different geometries, representing miscellaneous 
rotations. 

(a) 0i-fce[O,jr]. 

92 e [0, k] : 

( <P = <P2 + %, 
Q+ = 1 = 01 - 02, 

[f =-9i-§ (2tt) 

(b) fe-flie[o,4 

92 e [0, tc] : 
C 4> = 92 + f, 

flf! = { = 02-01, 

[9 = -9i-f (2n) 

(c) 0i + 02 e(0,*]. 

f * = 92+f, 

i2+ = ^ = 0, + 2 , (13) 
[f / =-9i + f (2jt). 

(d) 0! + 2 e [k,2k]. 

{ * = 92-f, 

I2 2 = ^ = 2a: - 0i - 2 , (14) 
l«P = -9i-f (2k). 

The function 12 j + = {I2 1 + 1 ,I2 1 ^} is matched for the case when Q\ > 2 , while the function 
= {£2^,£2^ 2 } ls matched for the case when 02 > 0i. The functions £2^ describe rota- 
tions for given 92 G [0, k] ; according to 12 j^, rotations are realized for (fa 6 (k, 2k] . The func- 
tion 12* = {^2^ , &2 } defines another possible rotation for the given angles 0i + 02 € (0, k] 
or 0] + 2 € [k,2k]. Note, if 92 £ [3^,2re], then for the rotation £2£, the angle <P > 2k. 
On the other hand, rotation over the angle 2k geometrically is equivalent to initial state. 
Thus, one can choose whether <P > 2k or < <£> — 2k < 2k. The same idea is valid and for 
the rotation Q2 ■ Finally, when 0i = 2 = 0, the angle = and the angles <£>, f acquire 



92 € (k,2k\ 



{ <*> = 92-f, 

£2+ = 1 = 0! - 02, 

I «P= -91 + f (2k) . 



(ID 



92 e (k,2k\ 



* = 92-f, 

= 02-0!, 

f = -9i + f (2») ■ 



(12) 
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any values in [0, 2k] . Then the rotation matrix R(Q) = R z (<£ + f) . Hence, in this case the 
full rotation is made by the angle <P + f around the z-axis. In other words, if Q\ = 8 2 = 0, 
rotation by the Euler angles Q is not singularly defined. Further it will be assumed that 6\ 
and 02 are not equal to zero at the same time (that is why the range of 0i + 02 for Q 2 is 
open from the left). 



4 Spherical functions 

In the previous section the mapping from S 2 x S 2 to SO (3) has been defined. It was demon- 
strated that possible rotations in R 3 from K\ to K 2 could be realized by the rotation an- 
gles £2f = {Q^,£2~} = {Q^ l ,Q^ 2 ,Q n ,Q^ 2 } and Q 2 = { Q.^ , Q. 2 } . If substituting these 
functions in Eq. (O, we would obtain the following spherical function (for alternative ex- 
pressions, see Appendix iBb 

(n,n^a,j3,7,5|x I ^2)J 9 ,=i a?+V (-l) 2(nt+ ' ,V )j3 9 '- 3 a(A; ) ?,?') (r(«*-«>0 

x {cos \\ (0! - 702)]} 2 *E^ (W) {tan[l (0! - yQ2)\? p + 4 ~ q , (15) 
p 

where Are Z+,Q+ and Q+ = {m+ 1/2; m 6 Z+}; the indices q,q' e [-k,+k\. Let us define 
particular cases of Eq. dl5b in the following way 

(0,n';+,±,+,-\xi,x 2 ) k qq i = ± %q q >(xuh), (16) 

(0,n';-,±,+,+\x h x 2 ) k qq , = ± tf q ,(x l ,x 2 ) = {-1)"'-* ± ^(x l ,x 2 ), (17) 

(0,n';+,+,-,+\x 1 ,x 2 ) k qq ,= + $ qf {x 1 ,x 2 ), (18) 

--,-\x u x 2 ) k w = -^(x u x 2 ) = (-l)^ + ^(xux 2 ), (19) 

and the matrices 

i2±: V^l.^Gi^^i,^),^^!,^)}. (20) 

It is seen, the spherical functions ± f] qq i {x\,x 2 ) and ± ^ k q i (xi,x 2 ) are the generalized spheri- 
cal functions D k qq , (i2 1 ± ) and D qq , (C2 2 ) parametrized by the coordinates of S 2 x S 2 , respec- 
tively. Particularly, the matrix + T\ k {x\,x 2 ) represents rotation in R 3 from K\ to K 2 , when 
01 > 02 and it is associated to the map £lf , while the matrix ~r\ k (jci, ^2) describes rotation 
from K\ to K2 when 02 > 0i and it is associated to the map Q p. The matrices + £* (x\,x 2 ) 
and ~ £, k (x\ , X2) are related to the maps Q 2 and Q 2 ■ 

In accordance with Eqs. (16M19I, the spherical functions ± E,, ± $ and ± £ are connected 
to each other as follows 

+ fik _ -Sk -rk _ +rk W 

V qq' ~ ^qq" V/ ~ ^ > 

It directly follows from Eq. <2U . that 
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+ T lqq' ~ ^qq' ~ ^qq' 6 { + < ^V' "=?«'}> ~ +< '-?-<?" ^2) 

where the column vectors of ± ^ k are orthonormal, i.e., 

I + & ^V 9 " = <V- (23) 
1 

The latter condition is, of course, valid and for the rest of spherical functions. 

It is noticeable, the matrices r\ k and ± ^ are the unitary irreducible matrix representa- 
tions of SO (3) (for k 6 Z + ) or of SU (2) (for k 6 Q + ), parametrized by the coordinates of 
S 2 x S 2 . Hence, the irreducible tensor operators T k transform among themselves as follows 
(to compare, see Eq. (0) 

A?: T k {K 1 )=Y f T ) k qq ,(x l ,x 2 )T k q ,(K l ), 

Of : T k (K 2 ) = £ ±QV M) m. (24) 

The transformation formula for the maps Q,f is restricted by the condition Q\ ^ 62 ■ In a 
contrary case, only the maps 12^ are valid. Consequently, the reduction formulas for the 
spherical functions 

{<'(•,•), %V (•>•)} ^ 

are these 

where in the brackets (.,.) the spherical coordinates of S 2 x S 2 are given. The Clebsch- 
Gordan coefficients of SU (2) are none zero only when q = q\ + qi and 4 = 4\ + ?2- The 
summation is performed over k = k\ — k 2 | , | k\ — ki \ + 1 , . . . , k\ + k 2 . 



k\ k% k 




k\ ki k 


q\ qi q 




q\ q' 2 4 



(26) 



Example. Suppose, x\ = (f ,f ) and x 2 = (j ,K). Possible rotations are realized then in 
accordance with X2j~[ = ? f >x)> ^2 = (iTi § > f) - Transformation formulas in Eq. < f24b 
are valid for the spherical functions ~l; k , ( f 1 f , § i ft) and + , ( ? , | , ^ , n) which coequal 
to D qq> (f . f . x) and (f . I • f ) 9 Suppose * = § and 9 = -±, = |. Then 

-«f.3 (?,f,f,*) =^.3 (f ,f ,? ) = (13-3^3) , 

5 5 5 

+ /-3 /7t n n _\ _ n 3 (M 2. E\ — 

S_ 1 3 (6' 4' 3' ;I 7 _ _A 3 I 2 ' 2' 4j _ 4 ■ 
3 3 3 3 

One may notice, according to expressions in Eq. ( 12 1 b . it is possible to find out how all other 
spherical functions are related to the calculated functions above. 



7 



5 The integral of spherical functions 

The choice of geometries £2 j , £2 2 i s convenient for other applications of obtained spheri- 
cal functions ± T , . This is because the parameters ft , ft and (p\ , 92 are separated. Conse- 
quently, the functions which depend on these parameters can be integrated separately, i.e., 
the integral on S 2 is naturally separated into the integrals over 9, and over ft By definition, 
the functions ^T k , (x\ ,$.2) are determined in different areas L 2 (£2) C S 2 , which restrict the 
existence of integrity. Therefore, we construct the integral 

■TtffalY) =1^2 {n,n';a,p,y,5\x h x 2 ) k qir (27) 

Suppose, 7= +1, i.e., the map Qf 1 is realized. Let the areas (or paths of integration) on 
S 2 to be 

L 2 (£2+) ={<p2e [0, n] ; ft e [0, ft]} , L 2 (£2+) = {(p 2 e [w,2w] ; ft e [0, ft]} , 
L 2 (i2 n ) ={<P2e [o,7t}-e 2 e [d u n]} , l 2 (n n ) = {92^ Mrc];fte [ft.rc]}. ^ ; 

Here £2^ and £2^ 2 mark the existence of integrable spherical functions ± ^ qq r {xi,X2) in the 
corresponding areas L 2 (i2 1 ± 1 ) and L 2 (£2^). The parameters a, /3, 7, 8, n are determined 
then according to the values presented in Table [Tj and the equations in Eq. d 10b . It is clear, 
the spherical function + £, is integrable in L 2 (£2^) with n 1 6 {1,2} and in L 2 (£2\ 2 ) with 
n' e {0, 1}; the function ~ E, is integrable in L 2 (£2^) with n' 6 {1,2} and in L 2 (fl^) with 
n' e {0, 1}. But -E, k , = (-iy-? , (see Eq. ©). This implies 

=/ i2K) ^ + ^^+J L2{Qn f 2 

= K q , (91 ) i^-i ^(— l) 9 ' + lj a{k,q, q 1 ) e~** 

x £> p (ft, ?1 </) (ft ;+;0, ft) + P I k qq , (ft ; +; ft , *)) , (29) 

f (-!)«', <Pie[0,f], 
V(<Pi) = < (-l) v , <Pie(f,f], (30) 
[(-1)3?', ^ e (3£,2»]. 

The definition of p fi, is given by the formula 

p I k qq , {d 1 ;r,a,b) = | 6 dft sin ft {cos [1 (ft - 7ft)]} 2 * {tan [> (ft - yd 2 )]} 2p+q ' ~ q . (31) 

For 7= +1 the integration is performed making the change of integrand z = tan [(ft — ft) /2], 
After some ordinary trigonometric manipulations it acquires the form 



D l k qq , (ft ; +; a, Z?) = 2 {a, b) cos ft + (/f+ (a, b) - 1% + (a, ft)) sin ft } , (32) 



ir(a,b)=lP + tan 



di-b 



7f+ tan 



ft -a 



(33) 
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with s = 0, 1 , 2 and /' (z) defined by 



/• Jlp+q 1 — q+s 2p+q' -q+s+\ 

7f+(z) = /dz^— 

JR (1+; 



(l +z 2)*+ 2 2/7 + ^-g + j+l 



,, / 2p + g , - g + ^+l 2p+V-g-M + 3 2 N 
2^1 ^ ,k + 2; ;-zr + const, (34) 



where 2F1 denotes Gauss hypergeometric function. 

When 0i = or Q\ = it, the function if (0, it) depends on infinite variables (z = ±00). 
For this reason bellow the boundary values of p 7* . (0;+;O, it) are presented. 

1 . = 0. Then p I k , (0; +;0, it) = 4/f + (-°°) (since /f + (0) = 0), where 

FW = t^B(k + l-p+^ +P+ ^) (35) 



with 5 being Beta function. 
2. = it. Then 

,4, (»;+;0,») = (-l)^'4 ? ,(0;+;0,^). (36) 
For example, it directly follows from Eqs. Pil l. Q6b . that 

h = J'dd sin 2k+l cos? 6 = [1 + (-1) 7 ] 7f (-) , (37) 

jfH = Ifl(*+l,I±lY (38) 

where the left hand side of I2 is the same integral discussed by Pinchon et al. [12| (Eqs. 
(19)-(20)), when developing rotation matrices for real spherical harmonics. 
When k 6 Z + , the integral of the spherical function r\ acquires the form 

S% (*;+) = Sqoltr* ((-l) 9 ' + l) a (k,0,q) ^ W<Pl 

xYfipfcM) rfyfa+'An). (39) 
p 

Hence, for integer k, the integral of (0,ra'; 0£,/3,+, 8\Jci, X2)g„i is none zero only if q = and 
q 1 is even. Particularly, {x\ ; +) = Alt. 

Finally, the integration of (n,n';a,fi,—,5\x\,X2) k qq i must be proceeded. For the maps 
Q.2 there is no difference which value of 0i and 62 is greater or less (or equal). It has been 
demonstrated earlier, that for Q.^ the condition O<0i + 02<?r(oc=j3 = 5 = +,n = O) 
must be satisfied, while for Qj the restriction is: it < Q\ + 02 < lit (a = /3 = 8 = — , n = 1). 
But when integrating ± £ k , (xi,x~2) over 02, the angle 02 acquires all values in [0, it]. This 

qq 

means the map 12^ is valid only if 0i = and the map Q.^ ls realized only for d\ = it. 
Hence, the integration of ± £ k , Ut\,x%) can not be correctly performed for any values of 

qq 

01, except for 0i = or 0i = it. In other words, the spherical functions ± £, k q i (^1^2), which 
represent rotations Sit , are not integrable on S 2 , in general. For this reason we conclude, that 
the most preferable spherical functions (at least for integration) are r\ k , {x\,x,2), i.e., those 
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which represent geometries Q^ 1 . Note, the angles (0i,(pi) and (62,92) are fully determined 
on S 2 for the maps Q,^ 1 . This implies that for any values of 6\, <pi, 02, <P2 there will always 
exist at least one rotation from the set Q, * = { Q, ^ , Q, ^ } . Consequently, the loss of geometry 
12* does not imply the loss of generality by means of the existence of at least one spherical 
function from the set { + E, , ~~ £ }. 



6 RCGC technique 

Since the irreducible matrix representations tj^ (#1,^2) depend on known coordinates x\ 
and Jt2, it is worth to exploit them in the study of tensor products of irreducible tensor 
operators T (or basis functions <j>) directly, but not formally, as in most cases, when tensor 
operators are transformed under representations (Q). Moreover, most of the physical 
operators 7]* basically studied in atomic spectroscopy, are expressed in terms of D and their 
various combinations. These are, for example, the spherical operators C*(i) = i k D^ Q (Q), the 
spherical harmonics Y^(x) = \J (2k + 1) /AnCz{£). Expressions over D of other operators, 
such as spin operator S 1 , angular momentum operator L 1 , can be found, for instance, in 
[8]. All these mentioned operators also transform among themselves according to Eq. ( 124b . 
Hence, it is natural for such operators to write T(x) instead of T(K), which is a general case. 
We say the operator T(x) acts on x coordinate. Below the latter notation will be used. 



6. 1 Transformation coefficients 

Reduction formula for tensor product reads (see Eq. J24H 

T% (*i)T% (h) =K {*i)<%*g n (iiM, (40) 
A 

In this paper the so-called for simplicity rotated Clebsch-Gordan coefficient of the first type 
or simply RCGC I is defined by 



m\ m'2 m 



4 1 if m ^,^2)=E^ ffl ,(x 1 ,X2) 

m 'l 

From Eq. J40t we obtain the next expression 

rX a \ .^AjAtA'A / a ^ \ ^pX' . 



(41) 



f * (* lt *z) = E (xi,x 2 ) t, (i! ) , (42) 

I'm' 



where rotated Clebsch-Gordan coefficient of the second type or simply RCGC II is 



CAiAtA'A /a a \ \ ' AiAtA' /a a \ 

! (XI, X7 = ) C ,\X\,Xl) 



m 1 w?2 



A] A2 A 
mi 7W2 w 



(43) 



Irreducible tensor operators T£ (x.\ ) and T£ (£i,X2) are delivered by applying reduction rules 
for the Kronecker product Ai x X2 — > A. However, T£ acts on x\, while 7^ acts on x\, x%. 
For example, if T represents normalized spherical harmonic C, then coupled tensor prod- 
uct t£ (fi) = {h+h-^c^ (x x ) ^ 0] ■ Contrarily, tensor operator f£ (x u x 2 ) can not be 

reduced into one due to different spaces, in which T m \ (xf) (i — 1,2) act on. 
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Another useful circumstance for RCGC application is based on a possibility to reduce 
these coefficients. In accordance with Eq. ( 126b . we directly attain 



(.,-)■ 



A 2 



"A, 


A 2 


X ' 




~\ 


h 


X 


mi 


m' 2 


m 




fh\ 


m' 2 


m 



X2 X2 A2 




X2 X2 A2 


mi ffi2 M2 




m' 2 m' 2 M' 2 



(44) 



Combining Eqs. d43b- (144b . we may also derive reduction formula for RCGC II. It is obvious, 
demonstrated reduction procedure (see Eq. ( 142b ) can be extended and for many-electron 
wave functions. 

Suppose, for instance, one needs to calculate matrix element of 7r (i) on the basis of 
eigenfunctions Y,', (x), described on SO(3)/SO(2). The application of RCGC technique di- 
rectly indicates, that in present case 



[l\\T k \\l\ = i>- l [(2l+l)/(2l + l)]V 2 T$(0) 



I k I 
000 



(45) 



The proof is produced in Appendix|C] Here it is assumed 6 = (0, 0). Reduced matrix element 
[l\\T k \\l\ is obtained from Wigner-Eckart theorem (lm\T k \!m) = (-l) 2k [l\\T k \\I\ [| * /„j . 

Particularly, if T k = C k , then operator Tq (0) = i*, and obtained formula is in fully con- 
sistence with Eq. (2.52) in | 2|. Thus, Eq. ( 145b is a generalization of this special case. Going 
on this route, one can derive matrix element expressions on the basis of more complex 
eigenfunctions. 

One can see from Eq. (142b . the specific feature of technique, based on coordinate trans- 
formations (or simply RCGC technique), is that tensor structure of operator and wave func- 
tions is preserved, and the resultant matrix element is calculated on the basis of transformed 
operators T (Jc), i.e., the calculation of multiple integrals transforms to the calculation of a 
single integral (over x). This is because each TV-electron matrix element 



dii / A.X2 ■ 
s 2 Js 2 



s 2 



dx/v 



ket 



(xi,X 2 ,...,*;v) 



XT Q (x U X2,...,X N )<P m , (X\,X2,...,X N ) 



is transformed to 



; / dx 2 / 

Js 2 J$ 



dx3 . 



libra A ' A ' 

&CN T]m 2 ll 2 {XI,X 2 ) • --rim N Nt l N {xi,X N ) 



Js 2 JS 2 
x (xi,x 2 ) r] k l K} (xi ,x 3 )... rj k ^ N (xi,%) 

yjeet jjcet yket 

x ^ (xuxi) n m , 3ll , (xi,x 3 )... r\J N ^ (i, ,%) , 

where <t> and T indicate coupled tensor products of ty^ (Jc\) and T k > (x\), respectively. The 

application of Eq. (f26]l for all TJ m |^ T\q i%% (xi,x%), TjJ^, (ii,^), £ =2,3,..., N 

implies, that initially determined 2./V-integral reduces to a double one 



^(*;+)...J^(<i;+)4# 
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Instead of that, for given iV-electron wave functions we produce N — 1 functions 5? (see Eq. 
d29l l) and a product of momenta coupling coefficients (CGC), which can be decomposed 
into 3n/-coefficients, if summing over projections. 



6.2 Example 

A simple application of RCGC technique can be demonstrated, for example, in a study of 
two electrons, located in some external field (of fixed nucleus, for instance). We refer to (9), 
where a two-electron wave function is presented by 



•Pm(n,r 2 )=^(r 1) r 2 ,0 I2 )Z)^ (Q), 



(46) 



with r, ■ = |r,|; / e Z + . In the same paper (Bhatia et al. \9\) it was determined, that Laplacian, 
involving dn (the angle between vectors r\ and r 2 ), does not affect the orbital angular 
momentum / > 0. Thus we may assume, that g is a radial function not going into a deeper 
analysis, since our aim is the angular part. Then f^' may be rewritten as follows 



fm (n ,r 2 ) = £^ 0i , r 2 , 6i 2 ) T] l mjl (*i ,x 2 ). 



(47) 



Suppose we want to calculate the Coulomb 1 jr\i matrix element. In this case the interaction 
in a tensor form reads T ° (r\,r2) = L*('"</'"> +1 ) {C k {xi) -C k (i 2 )), where r< = min(ri,r2) 
and r> = max(ri,r2). The scalar product (C (Jci)-C*(jc 2 )) is reduced in agreement with 
Eq. ( 142b . what leads to the following construction of (f^lr^ 1 f^') (taking into account only 
angular part) 

-l)" , - m 'E(- 1 )' 1 ^ ( r ^2, 0i2)^ +m _ m , (n,r 2 ,0 12 ) 



k k K 
00 

a y— even 

k k K 

m — m' Q + m' — m Q 

L 



-jX Q + m' — m Q + m' 



£^dxC*(*)jf e (x;+) 

I k L 
-m m — m' —m 

l it 

Q + m' — m — fi m — m'+jl Q 





L 


I' 


L 




—m 


m' 






(48) 



It is evident, this way of calculation is more efficient in comparison with direct integration 
of JJ S 2 diidi 2 D[ m (Q) (C k (ii) • C k (x 2 )) D^, , (£2), because otherwise one should change 
integrands, when mapping from x to Q . Consequently, that would lead to complex manipu- 
lations of trigonometric equations, given in Eq. ©. 



7 Summary 

In present study, we produced a parametrization of standard Wigner D-function on SU (2) by 
the coordinates (Jfi,Jc 2 ) of vector in the fixed and rotated coordinate systems. As a result, we 
found the set of spherical functions, conformed to D(Q ) in miscellaneous areas L 2 (Q ) c5 2 . 
We showed that offered parametrization of D(Q) provides an opportunity to reduce angular 
2,/V-integrals, which are of special interest in theoretical atomic spectroscopy, to a double 
one. Particularly, we demonstrated a new way of construction of irreducible tensor operator 
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matrix element, which plays the role of a generalization of previously obtained special cases 
(see, for example, Eq. d45ll). Another convenient usage of suggested RCGC technique, based 
on coordinate transformation, is applied to the calculation of matrix elements on the basis 
of functions, expressed in terms of standard D(Q) functions (for two-electron case, see Sec. 

1531 . 

Acknowledgements The authors are grateful to R. Kisielius for attentive revision of this paper and for 
various useful remarks. 



A Optimal values for f 

Here the solutions £2 = {<p + 7t/2,B,W) 6 RV.v; 6 S 2 of Eq. {6j will be studied in a more detail. Partial 
solutions of the systems (x2 , zi ) and (y>2 , zi ) for 9 , (p are 



zi z 2 + 02 |/i | Jzj-zi+y? 
Qo, <t-i = °~l arccos = — V 2ltn\ ,n\ GZ , (49) 



<Po 3 a t = 03 arccos ^ h27t«2, «2 e Z + , (50) 

1 — Zj 



<Pa 3 <T 4 = °3 arccos TTTT, r\ h 27T«2, «2 G ^ , (51) 



where o~, E { — 1 , + 1 }; sgn {y\ ) denotes the sign of y\ . Let us mark the common values of <pSj 4 and <Pa^a 4 by 
Wct 3 <7 4 (^b-?2;f)- Let us also denote 



2 2 (zf+3'i ) 

Obtained functions x and w ai'e substituted in Eq. (2)- The generalized spherical function D is rearranged to 
the following function 

{„ 2 Z* (jei,Jfe;¥)) ^(W^+^^^^LMW) 

t J 02 a 3 a 4 



xxS 2 £lt (.viA; l f')(l-x (l2 (i 1 ,.V2;f)) p+ *+V . (53) 

In order to look for optimal f values, we carry out variational procedure for the gauge f, applying it for all 
possible distributions of o~2, o~3, 04, i.e., (5/5f) n-,Z k qq * = 0. This implies 



: x. 



-^-'(l-x^)-'^^- 1 (p-Szt-^ ) =0. (54) 



Applying obvious fact, that A = B = OifA,BeRinA + iB = 0, we finally gain 

X02 =77,7,6 {0,1,^}, (55) 

It follows from equation x<j 2 = 7], that v'j 6 C. Thus, f does not belong to R. Equation -j^r = is rewritten 
in the form 

<5xo, _ <5x a2 Sy[ _ Sx a2 
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The solutions of gj = for do not belong to R. Thus, equation x\ = has to be solved. The latter is 

equivalent to Eq. Finally, when studying the third equation in Eq. i55\ , we would get some individual 
values of f , which depend onq.q' (for q = 0, present equation vanishes). Thus, the set of solutions would be 
the subset of solutions given by the first equation, independent of q,q'. 



B The alternative expressions of spherical functions 

In accordance with Eqs. {5J-j4), the product of coefficients a and b p can be rewritten as follows 



(k + q)\(k-q)\ 



(k + q')\(k-q>) 



*(k-4\( k+q 1 



p + q -q ) • 



(57) 



where the last two quantities on the right hand side of Eq. 1571 denote binomial coefficients. Further, let us 
mark the term 

A\ 4 (7) = a (k,q,q') {cos [■ (8, - Y^W^P (W) {tan [\ (9, - ye 2 )]} 2 " + "' '~" . (58) 



According to Eq. I57K A k , (y) may be revised by 



A k qcf (y) = v- q 



(k + q)\(k-q)\ 



(k + q')\(k-q')\ 



B-i)' 



p ) \p+q -q) (i+z) 



„(k-q'\( k+q' 



where z = tan 2 [i, (0i — 762)]- Performing the summation over p, we obtain 



„fk-q'\f k+q' \ z" + ^ 



P + 1 -1) (1+z) 



(59) 



= (-Vy-l (^J^Z^ {\+z) M 2F i {k + \+q,k+Y-c/-l + q-q>--z) (<?></) (60a) 

= ( K ] + q ' q y^(l+z) i+k 2F l {k+l+q',k+l-q-,\+q'-q-,-z) (q' > q) . (60b) 
Thus if: (a)A>* = A k ,, for q>d\ (b) A <k , =A k „foro' > q, then 

v ' qq' qq' ' 1 — ^ qq' qq' ' ^ — 



(y) = ^ [ MM 

l«' {n (q-q')\ [{k + q')\{k-q)\ 



{tan [\ (9, - r&2)]} M {cos [J (0! - yfe)]}" 



x 2 Fi (A+l + g,fc+l-« / ;l+ ? -?';-tan 2 [^(e 1 -7e 2 )]) . 



^(r)=A^(r). 



Substituting A in Eq. 115) we gather 



(n,n;a,B, 7,5*1, xi) , = < , , „ ' '*£ 
(n,B^a,j3,7,5|ii^ 2 ) > j=i a?+ ^(-l) 2 ^ + " V )jS«' _ V( ? ' R 2-«'« , 0A^(r) ) 
(n,n / ;a,P ) y,«|jei ) jfe)^=i( - a K«-4') e i(ft+fl!)(«-«') (n,n';a,p, y,8\ Xl ,xi) >k ■ 



(61) 
(62) 

(63) 
(64) 
(65) 
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C Reduced matrix element on the basis of spherical harmonics 

Here a proof of Eq. 1451 . exploiting RCGC technique, will be offered. The matrix element of 7 1 * (x) on the 
basis of arbitrary functions Xjf'J (x) and iff^ (x) is written by 



I k I 
m q m 



(lm\T q k \lm) = [l\\T k \\l] 
Direct adaptation of Eq. {24) points to 

(lm\T*\lm) = Jit I ^;^)^^)^^)^^^)^^^)^^'^)- 



(66) 



(67) 



Transformed functions \f/ and operator T depend on fixed coordinates x 1 . Consequently, they can be located 
in front of the integral. Recalling, that Tj' ; = (— l) m ~ m 7j' _ and exploiting reduction rules for the Kro- 
necker products / x & — > L, L x f— > L (see Eq. (26)), we attain rf^ u i (x' ,x). The integral of present spherical 
function is defined in Eq. j29h and in this case equals to •9"n M i{^',+) (see Eq. (39)), since L 6 Z+. Thus 

(lm\T k \lm)= £ (-l) m -"'V^(^)^(x')^oV(^;+)^(^) 



/ k L 




I k L~ 




L I L 




L I L 


—m q —m 




-m' q' M' 




—m m 




M' m' M' 



(68) 

The matrix element does not depend on x', thus one can choose any value. We select the most simple case 
x! = = (0,0). The application of orthogonality condition for Clebsh-Gordan coefficients, combining Eqs. 
(66), (68) , leads to equalities L = J, L = 0; thus .5^° (0;+) = An and 



\\T k \W 



4n 



21- 



T Eyj(6)^(6)v4(8) 



/ k i 



(69) 



It is noticeable, obtained expression, by applying technique of coordinate transformations, coincides with 
Eq. (41) in | 8 |, if reducing given Kronecker products / x k — > L', L' x /~ — > L". This would signify L' = I 
and L" = 0. At this step we turn to special case of eigenfunctions V'mt*) = ^m(^)' what implies equality 
^m(O) = 5„,oi' -\/(2' + I) /An. Hence Eq. (69) becomes equal to expression, presented in Eq. 145 1 . One should 
be also reminded, that in general, in Eq. 169 \ the arguments can be replaced by any values x. 
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